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Abstract 

We have performed measurements of the magnetic susceptibility of NaA^Os 
between 2 and 400 K. The high temperature part is typical of spin 1/2 chains 
with a nearest-neighbour antiferromagnetic exchange integral J of 529 K. 
We develop a model for the susceptibility of a system with finite chains to 
account for the low temperature part of the data, which cannot be fitted by 
a standard Curie- Weiss term. These results suggest that the next nearest- 
neighbour exchange integral J2 in CaV^g should be of the order of 500 
K because, like J in NaV20s, it corresponds to corner sharing VO5 square 
pyramids. 

PACS Nos : 75.10.jm 75.40.Cx 75.50.Ee 
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The report by Taniguchi et al |J of a spin gap behaviour in the quasi-two dimensional 
system CaA^Og has triggered an intensive theoretical activity aimed at understanding the 
origin of this gap [§H§|- The emerging picture is that there is no spin gap in the model with 
only exchange integrals J\ between nearest neighbours @,§,§J , and that there is a spin gap if 
a coupling constant to second neighbours J 2 is included as long as 0.2 < J2/J1 < 0.7 [[[],(|||. 
To check this theory, one needs information on the value of the exchange integrals. This 
information turns out to be difficult to extract from the suceptibility. The best calculation 
of the temperature dependence of the susceptibility of that model is a high temperature 
expansion due to Gelfand et al ||. Assuming J2/J1 = 1/2, they could reproduce the maxi- 
mum of the susceptibility around 100 K with J x ~ 200 K. The fit of the high temperature 
part is not satisfactory however, and the question of the value of the integrals is still pretty 
much open. 

CaX^Og is actually a member of a large family of vanadium oxides studied by Galy and 



coworkers in the mid seventies []10| , and a natural idea is to look at other members of the 
family to try to get information on the exchange integrals. The other 2D compounds that can 
be synthesized with Ca, CaX^Os and CaVaO^, lead to a similarly difficult problem because 
they involve both Ji and j \. This difficulty can be overcome by studying an other mixed 



is 



valence vanadium oxide, NaX^CV This compound, first synthesized by Hardy et al [JTT 
isostructural to CaX^Os. Note however that NaV20s crystallizes with the orthorhombic non 
centro-symmetric space group P2mn, while CaX^Os crystallizes with the centro-symmetric 
space group Pmmn. Now, NaX^Os contains Na + instead of Ca 2+ , and half the vanadium 
have to be in the oxydation state V 5+ (formally one has NaV 5+ V 4+ Os). These ions do not 
carry a spin, while remaining V 4+ carry a spin 1/2 and form a set of well separated chains 
of corner sharing VO5 square pyramids (SP) (see Fig. 1). The magnetic properties should 
thus be well described by the one-dimensional spin 1/2 Heisenberg model: 

H = JJ2S t .S l+ i (1) 

i 

Note that the exchange integral J between corner sharing VO5 SP is equivalent to the 



next-nearest neighbour exchange integral J2 of CaX^Og. 

In this Letter, we present measurements of the magnetic susceptibility of NaX^Os from 
2 to 400 K. The compound was prepared starting from a stoichiometric mixture of NaVC>3 
(Merck, min 99%), V2O3 (obtained by hydrogen reduction of V2O5 at 800 °C) and V2O5 
(Aldrich Chem. Co. ,99, 9%). The mixture was ground intimately, sealed in an evacuated 
quartz tube and then heated for 10 hours at 615 °C. The X-ray diffraction pattern of the 
resulting dark powder indicated the formation of the pure phase NaX^Os. A refinement 
of the structure of NaX^Os was proposed by Carpy and Galy ||12|| . It is orthorhombic 



and consists, as shown on the perspective view in fig. la), of two dimensionnal layers of 
VO5 SP with the Na atoms between the layers. It is worth mentionning in this structure the 
ordering of the V 4+ and V 5+ atoms in the layers with formation of rows (Fig. lb)). Magnetic 
susceptibility measurements were performed using a SQUID susceptometer. The magnetic 
field intensity was lkG. The molar susceptibilities were corrected for diamagnetism by using 
Pascal's constants. 

The raw data are presented in Fig. 2. They agree with the early measurements between 



80 and 600 K by Carpy et al [I3|. Above 100 K, the suceptibility is consistent with that 



of a spin 1/2 chain flT^JTo]] . In that temperature range, the best available estimate of the 



susceptibility due to Eggert et al |L5| is actually indistinguishable from the Bonner-Fisher 



result fli~4| |, so it does not matter which theory we use to fit the data. There is a maximum 



at 350 K which implies an exchange integral J ~ 529 K. 

Below that temperature, there is no evidence of a phase transition or of three dimensional 
ordering, but, as usual, there is an increase of the susceptibility at low temperature due to 
some kind of defects. The standard procedure is to describe these defects by a Curie- Weiss 
term x° W (T) = g 2 p 2 3 S(S + l)/3k B (T — 9), so that the spin part of the susceptibility reads: 

X tot (T) = (1 - p) X ™{T) + PX CW {T) + x VV (2) 

X VV is the temperature independent Van Vleck paramagnetic susceptibility, p is the concen- 
tration of impurities, and x°°{T) is the susceptibility of the infinite chain. At low tempera- 



ture, the difference between the Bonner-Fisher estimate and the recent results of Eggert et 
al [[15] is not negligible ||16|| , and we have used the results of Eggert et al for x°°(T). It turns 



out that the low temperature part of the susceptibility cannot be fitted satisfactorily along 
these lines. The best fit one can get using Eq. (^|) is depicted as a dashed line in Fig. (2). 
It considerably overestimates the actual susceptibility around the minimum at 70 K. 

The main problem is that the amount of impurities one needs to interpret the low tem- 
perature susceptibility gives a much too large contribution at higher temperatures. In other 
words, the susceptibility behaves as if the impurities were slowly disappearing when the 
temperature increases. While this clearly cannot be reconciled with extrinsic impurities, 
such a behaviour actually makes sense if the impurity contribution comes from finite chains 
with an odd number of sites. The idea is the following: Roughly speaking, a finite-length 
chain with N spins behaves like an infinite one at temperatures larger than the finite-size 
gap, and like a finite one below that temperature. Now, the finite size gap is of order J/N. 
So if we have a distribution of finite chains with different lengths, they will progressively 
dissappear from the impurity term to contribute to x°°(T) as the temperature is increased. 

To be more quantitative, we need to know the distribution of length of the finite chains. 
If we make the reasonable assumption that the finite chains are due to a random distribution 
of point defects, then elementary statistical mechanics shows that the distribution of length 
is of the form 

P(N) = EiEt^O (3 ) 

where the partition function Z = (1 — exp(— A)) -1 while the Lagrange parameter A is fixed 
by the average number of sites Nq of the chains according to A = log(l + 1/Nq) ~ I/Nq if 
JVo > 1. The average number of sites of the chains is itself related to the concentration p 
of defects by N — 1/p. Then the total susceptibility per spin is simply given by xCO = 
(1/iVo) J2n P(N)xn(T) where Xn(T) is the susceptibility of a chain of length N. 

Numerical calculations of xn{T) are indeed possible |TJ|], but even with today's numerical 



facilities, accurate results would be limited to relatively short chains. However, a good 



estimate can be obtained in the following way. At low enough temperature, a finite chain 
with an odd number of sites behaves like a spin 1/2 impurity according to Curie's law 
X C '"" e (T) = g 2 fi 2 B S(S + 1)/3/cbT, while the susceptibility of a chain with an even number of 
sites vanishes exponentially. Averaging the susceptibility according to Xn = (Xn +Xtv+i)/2, 
and neglecting the exponentially small contribution of the even chain, the susceptibility of 
finite chains is given by Xn(T) = x Curte {T)/2 at low temperature and Xn(T) = Nx°°(T) 
at high temperature. So, if we define the cross-over temperature T C (N) as the temperature 
where these two expressions are equal, we can approximate Xn(T) by x rie (T)/2 if T < 
T C (N) and Nx°°(T) if T > T C (N). Equivalently, one can define a length scale L(T) by 
fcsT c (L(T)/a) = J, where a is the lattice parameter. Then, for a given temperature T, 
Xn(T) will be given by iVx°°(T) if iV > L{T)/a and by x Cune (T)/2 if N < L{T)/a. Solving 
k B T c (L(T)/a) = J for L(T) yields L{T)/a = x Cune (T)/2x°°( y T). This expression can be 
cast in a more transparent form if one writes x°°(T) = ^j 3 ^ X~{T), where x(T) is the 
normalized susceptibility that equals 1/tt 2 at T = 0. In terms of this function, one has 
L(T)/a = (l/8x(T)) J//cbT. Note that 8%(T) is of order 1 for T < J. The cross-over 
temperature T C (N) ~ J/N is thus of the order of the finite size gap, as it should. 
So, the total susceptibility per site is given by 

( Curie (rn\ \ 
X / ] E P(N) + X°°(T) E NP(N)\ (4) 
N<L{T)/a N>L{T)/a J 

This equation is actually valid for any distribution of chain length P(N). Concentrating 
on the distribution of Eq. (§), the sums are readily performed: J2n<l(t)/ci P(N) = 1 — 
exp(-L(T)/aN ) and E N >L(T)/a N P(N) = N - (N + L(T)/a) exp(-L(T)/aN . Including 
a VanVleck contribution, our final result for the susceptibility reads 

X M {T) = l-*M-m/aN ) xCW{T) + + eM-L(T)/aN ) X °°(T) + X VV (5) 



with L(T)/a = (l/8x(T))J/k B T 0. We have replaced x CuHe {T) by x° W (T) to account 
for possible residual interactions between the finite chains. This can be seen as an extension 
of Eq. (H). The concentration p = 1/Nq is still an adjustable parameter. However, the 



coefficients in front of x CW (T) and x°°(T) are no longer simply p and 1 — p, but functions 
of temperature. Using Eq. ([5]), we have been able to obtain a much better fit of the raw 
experimental data. This fit is shown as a solid line on Fig. (2). The parameters are J = 529 
K, g = 2.043, x VV = U1.10- 6 ccm/mole, N = 35 and 9 = -1.26 K. This value of iVo 
corresponds to a concentration of defects p = 2.9 % which seems to be a reasonable number. 
The fit was realized with the help of a simplex non-linear least-squares fitting procedure, 
and the relative deviation defined as Ei(x° 6s - X calc ) 2 / 12(xt s ) 2 was equal to 3.1(T 5 . 

Let us come back to CaX^Og for a moment. Our data for the susceptibility of NaX^Os 
are consistent with the model of Eq. ([!]) with an exchange constant J = 529 K. This 
value should be contrasted with the values proposed for J 2 in the case of CaV 4 9 , which 
range from 50 to 100 K H0,f|. While one cannot exclude some dependence of the exchange 
integrals on the overall chemical environment, they should essentially depend on the local 
geometry, which is the same in NaX^Os and CaX^Og for corner sharing SP. Besides, the 
present determination of J is quite unambiguous because this is the only exchange integral 
involved. So we are led to the conclusion that previous estimates of J 2 in CaX^Og cannot 
be correct. If we put all the information we have at the moment, a consistent picture of 
exchange integrals can still be obtained. The presence of a gap of 107 K in CaV/iOg, i.e. 
much smaller than J 2 , can still be explained by the J± — J 2 model on the depleted lattice if 
J 2 /Ji is not too far from either of the critical values 0.2 and 0.7 where the gap disappears. 
Now, according to a recent work of Kontani et al |18[ , the observation by neutron scattering 



of stripe order in the compound C&V3O7 [|l£| implies that J 2 /J± cannot be too small. At a 
quantitative level, the bound given by the modified spin-wave theory J 2 /J\ > 0.6932 cannot 
be taken too seriously, but a ratio J 2 / J\ close to 0.2 can be excluded because CaN zOj should 
then exhibit Neel order. So we think that coupling constants Ji ~ 700 K and J 2 ~ 500 
K are the best candidates so far to describe this type of vanadium oxides. A very useful 
check will be to see if such values are compatible with the temperature dependences of the 
susceptibility reported for CaV 3O7 [^0J and CaX^Og JT|. 

In conclusion, we have developped a theory to describe the low temperature susceptibility 
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of a one- dimensional spin 1/2 antiferromagnetic system with finite chains. This theory 
leads to a very nice fit of the data we have obtained for NaX^Os. More generally, this 
theory should provide a much more accurate way of substracting the contribution of finite 
chains from the experimental data than just adding a Curie- Weiss term. This might help 
identify in other compounds the remarkable temperature dependence predicted by Eggert 
et al for the infinite spin 1/2 antiferromagnetic chain and already observed in S^CuOs 
according to Eggert's reinterpretation of the experimental data fl6| . It would be quite 
interesting to see how accurate the present theory is by performing Monte Carlo calculations 
of the temperature dependence of the susceptibility of finite chains. It would also be very 
interesting to test it on systems where the concentration of non magnetic defects that break 
chains can be controlled. Work along these lines on various vanadium based oxides is in 
progress. 

We acknowledge useful discussions with M. Albrecht, J. P. Daudey, M. Luchini, D. 
Poilblanc, J. -P. Renard and J. M. Savariault. We are especially grateful to J. Galy for 
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FIGURES 

FIG. 1. Structure of NaA^Os. a) Perspective view in the [010] direction. The square pyramids 
occupied by V 4+ are indicated by an arrow, b) Schematic representation of the (V20s)n rows 
along the Oy axis. 

FIG. 2. Thermal variation of the molar magnetic susceptibility of NaV20s. Diamonds: Ex- 
perimental data; Broken line: Fit using Eq. (2); Solid line: Fit using the present theory (Eq. 5). 
Insert: Enlargement of the low temperature region. 
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